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The  present  method  is  based  on  the  extension  of  the  Vortex  Lattice  concepts  also 
to  low  aspect  ratio  wings  and  the  rolled-up  vortex  which  is  established  over  as 
well  as  behind  the  wings.  The  computer  program,  which  was  developed,  is  capable 
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of  reproducing  the  results  of  the  Vortex  Lattice  Method  (vm)  for  linear  aero- 
dynamic coefficients  for  high  aspect  ratio  wings.  The  wake  effect  is  included 
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measured  data. 

Numerical  problems  and  limitations  of  the  present  method  are 
discussed.  The  finite  thickness  of  the  vortex  wake  is  being  studied 
by  using  a two-dimensional  method  and  assuming  that  the  wake  cross 
section  contains  vorticity  in  an  otherwise  irrotational  field.  The 
wake  development  under  its  own  induced  field  is  calculated,  ignoring 
viscous  dissipation.  A finite  spiral  wake  structure  is  observed  to 
develop. 

The  research  program  included  also  some  investigations  of 
viscous-inviscid  interactions,  particularly  problems  which  are 
important  for  wing-body  aerodynamics.  Calculation  programs  based 
on  various  iterative  processes  between  the  solutions  of  the  viscous 
flow  regions  and  the  essentially  inviscid  regions  have  been  studied. 
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I.  INTRODUCTION 

The  present  report  summarizes  the  research  program  conducted 
under  Grant  AF  OSR  71-2145  during  the  period  1971-1975.  Results 
of  this  research  program  are  presented  in  the  reports  and  papers 
listed  in  References  1 to  10. 

The  research  program  was  aimed  at  the  investigation  of  the 
rolling  up  of  vortices  over  and  behind  wings  of  various  planforms, 
and  the  simultaneous  evaluation  of  the  spanwise  and  chordwise  load- 
ings. This  unified  analysis  is  the  main  feature  of  the  calculation 
method  presently  developed.  On  one  hand  this  method  is  based  on  the 
representation  of  the  vortex  distribution  over  the  wing  planforms  by 
a vortex  lattice  similar  in  some  respects  to  the  method  developed 
by  Falkner  (Refs.  11,  12,  13);  on  the  other  hand,  the  effects  of  the 
vortices  from  each  cell  which  are  shed  into  the  flow  field  are  also 
considered,  in  the  present  investigation.  These  vortices  are  first 
allowed  to  leave  the  wing  planform  from  its  trailing  edges,  as  in  the 
investigations  of  Westwater  [14] , Hacket  and  Evans  [15]  and  Hancock 
and  Butter  [16] , Rom  and  Zorea  [2] . In  further  steps  the  three- 
dimensional  geometry,  the  shedding  of  vortices  from  the  side  and 
leading  edges  and  from  every  cell  in  the  planform  and  the  combined 
calculation  of  the  vortex  wake  and  the  wing  lift  are  developed. 

Since  the  present  work  makes  no  assumptions  concerning  the  trail- 
ing vortex  behaviour  once  it  has  left  the  wing  - assumption  which 


would  limit  the  calculation  to  a particular  range  of  aspect  ratios  - 


and  since  separation  of  the  vortices  is  permitted  at  all  edges  or  over  the 
whole  wing,  the  present  approach  will  cover  all  the  range  of  aspect  ratios 
in  an  unified  manner  provided  it  can  be  shown  that  separation  is  automatical- 
ly suppressed  when  it  is  inappropriate,  as  for  example,  in  the  case  of 
the  leading-edge  and  surface  of  a high-aspect-ratio  wing.  The  calculations 
of  Reference  5 show  this  to  be  effectively  done. 

Various  relevant  investigations  are  presented  in  Refs.  17  to  23,  and 

i 

in  the  publications  under  the  present  Grant  in  References  2,  4,  5.  Another 
important  aspect  of  the  vortex  wake  investigation  is  the  process  of  rolling 

I 

up  of  the  vortex  sheet.  This  process  was  first  investigated  by  Kaden  [24] 

and  by  Spreiter  and  Sacks  [25]  . The  effects  on  the  process  of  rolling  up 

I 

of  the  finite  thickness  of  the  wing  and  its  wake  was  investigated  as  part 
of  our  program.  The  effects  of  wing  thickness  on  the  flow  field  are 
presented  in  Refs.  1 and  7 and  an  investigation  of  the  initial  rolling 
up  of  a thick  vortex  layer  is  presented  in  Ref.  10. 

An  additional  topic,  investigated  by  both  lifting-line  theory  and 

vortex-lattice  calculations,  is  the  manner  of  the  effect  of  a number  of 

vortex  suppressing  devices,  which  have  been  tried,  on  the  overall  aero- 

dynamic  properties  of  the  wings  to  which  they  are  attached  [6] , [8] . The 

results  for  a Boeing  747  wing  fitted  with  eighteen  pairs  of  triangular 

vanes  for  vortex  suppression,  which  were  calculated,  as  reported  in 

Reference  8,  are  compared  there  with  experimental  results  given  in  Ref. 

33.  Good  agreement  is  obtained  for  the  reduction  of  C at  cruise  con- 
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dition  but  AC^  is  much  larger  than  that  given  by  experiment.  Most  of 
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the  CD  increment  found  in  the  calculation  is  due  to  the  vane  drag  and  not 
to  the  change  of  wing  induced  drag. 

An  important  result  of  the  new  method  of  References  2,  4 and  5 is 
the  evaluation  of  the  non-linear  aerodynamic  coefficients.  There  are  a 
number  of  methods  proposed  for  the  evaluation  of  the  non-linear  lift  of 
wings.  Some  of  the  investigations  that  form  the  basis  for  methods  of 
calculation  of  non-linear  aerodynamic  lift  coefficient  were  presented  by 
Brown  and  Michael  [26],  Mangier  and  Smith  [27,  28,29],  Gersten  [30]  and 
Polhamus  [31] . The  method  developed  in  the  pesent  research  make  the 
evaluation  of  the  non-linear  characteristics  of  wings  possible  as  sum- 
marized in  Ref.  5. 

The  computer  programs  developed  in  References  2,  4 and  5 follow  the 
procedure  by  which  the  wing  is  first  divided  into  cells.  Then  a bound 
vortex  element  is  placed  at  the  "quarter  chord"  of  each  cell  and  a control 
point  is  placed  at  the  middle  of  its  "three  quarter  chord".  The  trailing 
vortices  can  be  taken  to  be  shed  away  from  the  edges  of  each  bound  vortex 
and  to  continue  in  the  wing  plane  to  the  trailing  edges.  In  this  case 
the  variation  of  the  aerodynamic  characteristics  v-  . j essentially  linear- 
ly with  the  angle  of  attack  [2] . On  the  other  hand  these  vortices  can  be 
taken  to  leave  the  wing  plane  from  only  the  side  and  trailing  edges  or 
from  each  cell.  When  the  trailing  vortices  are  allowed  to  leave  from  the 
side  and  trailing  edges  only  the  trajectoties  of  these  trailing  vortices 
are  determined  by  an  iterative  calculation  which  results  in  the  evaluation 
of  non-linear  aerodynamic  characteristics.  In  the  case  when  the  trailing 
vortices  are  allowed  to  leave  the  wing  plane  at  each  cell  at  an  arbitrary 


I 


4 - 


{ 

* 


I 


I 

f 


local  angle  (similar  to  Gersten's  model)  these  vortices  can  be  either 
straight  lines  ( as  in  Gersten's  model)  or  free  to  roll  up  above  and 
behind  the  wing  planform  non-linear  aerodynamic  wing  characteristics  are 
also  obtained  [2] . 

The  numerical  problems  involved  in  approximating  the  circulation  dis- 
tribution by  discrete  ideal  vortices  have  been  investigated.  The  numerical 
difficulties  associated  with  using  discrete  ideal  vortices  become  even  more 
difficult  when  the  vortices  are  shed  from  the  cells  in  the  wings'  plan- 
form  or  from  its  side  edges.  In  order  to  eliminate  some  of  these  difficul- 
ties it  was  proposed  in  Ref.  5 to  use  the  finite  core  vortex  model.  As  a 
first  approximation  the  core  size  is  related  to  the  local  vortex  strength 
by  the  method  proposed  by  Sprieter  and  Sacks  [25] . 

An  additional  subject  of  investigation  dealt  with  the  improving  the 
linearized  solution  for  the  flow  near  the  tips  and  wake  edges  of  a lifting 
wing  in  incompressible  flow,  mainly  with  a view  to  its  more  effective 
utilization  in  the  calculation  of  the  rolling  up  of  the  trailing  vortex 
wake  (Ref.  1) . The  linearized  solution  always  exhibits  singularities 
of  velocity  at  the  wing  leading  edge  and  the  wake  edge  due  to  the  re- 
placement of  these  actually  or  effectively  radiused  edges  by  sharp  ones  in 
the  linearisation  process.  The  present  treatment  takes  account  of  the 
actual  or,  in  the  case  of  the  wake,  effective  radius  of  the  edge  by  the 
method  of  matched  asymptotic  expansions  and  so  eliminate  this  unrealistic 
feature . 
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The  result  takes  the  form  of  an  integral  for  the  perturbation  potential 
which  is  uniformly  valid  to  first  approximation  in  all  parts  of  the  flow 
field  about  the  wing  and  wake.  The  calculation  by  this  method  of  the 
potential  function,  the  velocity  and  the  pressure  distributions  on 
ellipsoidal  wings  at  zero  incidence  is  presented  in  Ref.  7. 

Investigations  of  wing  wake  roll  up  have  assumed  the  wake  to  be  a 
vortex  sheet  of  zero  thickness.  We  have  initiated  a study  assuming  that 
the  wake  cross-section  has  a finite  thickness  and  some  plausible  shape. 

A two-dimensional  method,  analogous  to  that  of  Westwater,  is  developed, 
assuming  that  the  wake  cross-section  contains  vorticity  in  an  otherwise 
irrotational  field. 

The  wake  is  divided  into  triangular  elements  and  the  vorticity  in  these 
is  determined  by  assuming  a linear  transverse  velocity  profile  in  the  wake 
and  that  the  initial  unrolled  wake  moves  downwards  as  determined  by  the 
wing  spanwise  loading  through  ordinary  wing-wake  theory.  Euler  time  step 
integration  is  then  used  to  calculate  the  wake  development  under  its  own 
induced  velocity  field,  ignoring  viscous  dissipation.  A finite  spiral  wake 
structure  is  observed  to  develop.  Some  preliminary  results  are  presented 
in  Ref.  10. 

The  research  effort  started  under  Grant  AF  OSR  71-2145  included  also 
an  investigation  of  viscous-inviscid  interaction;  particularly  emphasizing 
areas  and  problems  which  are  important  for  wing  aerodynamics.  The  inter- 


active flow  analysis  is  based  on  dividing  the  flow  field  to  regions  where  the 
flow  is  essentially  inviscid  and  regions  where  the  flow  is  dominated  by 
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II.  THE  GENERALIZED  VORTEX-LATTICE  TECHNIQUE  - METHODS  OF  CALCULATION 

OF  THE  TRAJECTORIES  OF  THE  VORTICES  AND  THE  AERODYNAMIC  CHARACTERISTICS 

OF  THE  WING. 

The  methods  of  calculation  are  based  on  a unified  procedure  for  cal- 
culation of  vortices  distributed  on  the  wing  planform  and  the  vortices  shed 
from  the  wing  into  the  flow  over  and  behind  the  wing.  The  results  of  this 
calculation  give  on  one  hand  the  loading  distribution  on  the  wing  planform 
which  makes  possible  the  evaluation  of  the  aerodynamic  characteristics 
(such  as  lift,  induced  drag,  aerodynamic  moments  etc.)  and  on  the  other  hand 
this  calculation  is  used  for  the  determination  of  the  flow  field  over  the 
wing  and  the  trailing  vortex  flow  field  behind  the  wing,  so  that  the  wake 
geometry  may  be  iteratively  relaxed  into  the  proper  form. 

The  calculation  procedure  follows  the  following  steps  (shown  schema- 
tically in  Fig.  1) . 

2.1.  Geometrical  Division  of  the  Planform  into  Cells 

The  wing  planform  is  first  divided  into  cells  by  proper  spanwise  and 
chordwise  divisions.  This  division  is  similar  in  many  respects  to  that 
used  in  the  standard  vortex  lattice  method.  One  has  to  select  a reasonable 
number  of  divisions.  A small  number  of  divisions  (large  cells)  will  result 
in  insufficient  accuracy  while  an  excessive  number  of  divisions  will  cause 
numerical  problems  (as  discussed  in  (Ref .2)  ) . Each  cell  is  then  represented 
by  a separate  horseshow  vortex  with  the  bound  vortex  positioned  at  the  cell ’ s 


1/4  chord  line  and  the  trailing  vortices  are  either  embedded  in  the  wing 
planform  and  leave  the  wing  at  the  trailing  line  or  are  free  to' leave  the 
wing  planform  right  at  the  edges  of  the  bound  vortex.  It  is  also  possible 
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to  allow  only  the  vortices  from  the  leading  and  side  edges  to  be  free 
to  leave  into  the  flow  while  all  other  trailing  vortices  are  embedded 
in  the  wing  planform  and  leave  the  wing  at  its  trailing  edge. 

The  vortex  strengths  are  calculated  to  be  compatible  with  the  boundary- 
conditions  applied  at  control  points  positioned  at  the  middle  of  the  3/4 
chord  lines  of  each  cell.  For  rectangular  wings  the  number  of  chordwise 
and  spanwise  divisions  can  be  selected  as  required,  whereas  in  the  case 
of  delta  wings,  it  is  found  that  we  must  use  equal  number  of  spanwise  and 
chordwise  division  lines.  The  bound  vortex  lines  and  the  control  points 
are  positioned  as  shown  in  Fig.  2,  where  examples  of  the  divisions  of 
rectangular  and  delta  wings  are  shewn. 

2.2.  Selection  of  Appropriate  Initial  Vortex  Distribution 

The  initial  vortex  distribution  is  determined  by  the  physical  model 
used  for  the  flow  field.  These  models  can  be: 

(1)  All  vortices  embedded  in  the  wing  planform  and  leaving  the 
wing  at  the  trailing  edges  into  the  wake  and  rolling  up 
behind  the  wing. 

(2)  Trailing  vortices  leaving  the  1/4  chord  line  of  each  cell  on 
the  leading  and  side  edges  only  ( rectangular  planform)  or  from 
leading  edges  only  (delta  planforms) . The  other  vortices  are 
embedded  in  the  wing  planform  and  leave  at  the  traling  edges 
only.  The  rolling-up  starts  at  the  side  edges  or  leading  edges 


for  rectangular  and  delta  wii.js  respectively. 


vortex  pattern  is  calculated  step  by  step) . The  correct  combinations  were 
found  by  a series  of  numerical  experiments.  It  is  important  to  note  that 
due  to  the  causes  mentioned,  too  many  vortices  or  stations  can  lead  to  the 
deterioration  of  the  results  rather  than  convergence  to  a limiting  form. 

In  the  next  stage  of  the  work  a wing  with  an  elliptic  circulation 
distribution  is  dealt  with,  employing  a single  bound  vortex  on  the  wing 
quarter-chord  line  (see  (b)  in  Fig.  3) . The  trailing  vortices  are  con- 
tinued straight  from  this  line  to  the  trailing  edge  and,  from  the  trail- 


ing edge  downstream.  Their  rolling-up  is  taken,  as  a first  approximation, 
to  be  given  by  the  results  of  Section  (a) . The  corrections  necessary  to 


this  wake  shape  to  conform  with  the  lifting-line  vortex  pattern  are  then 
calculated  by  a series  of  iterations,  assuming  the  elliptic  spanwise 
circulation  distribution  unchanged  throughout.  The  experience  gained  in 
stage  (a)  is  utilized  to  select  the  number  of  vortices  etc.  At  a later 
stage  the  finite-core  concept  is  incorporated  here  too. 

A practical  point  worth  mentioning  is  as  follows:  The  rolling-up 

calculations  only  carried  as  far  as  ten  semi-spans  downstream  of  the 
trailing  edge.  However,  a correction  for  the  influence  of  the  portion  of 
the  wake  downstream  of  this  station  was  incorporated.  At  the  ninth  and 
tenth  stations  the  center  of  gravity  of  the  vortex-wake  cross-section 
(in  the  half-span  domain)  was  found  and  a line  through  these  two  points 
was  taken  to  define  the  position  and  direction  of  a single  replacement 
trailing  vortex  from  each  half  wing,  starting  from  the  tenth  station 
extending  to  infinity  and  having  the  appropriate  ultimate  strength  (i.e. 


equal  to  the  value  of  the  mid-span  circulation  - similar  to  that  used 
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in  [15]).  The  induced  velocity  due  to  this  pair  of  semi-infinite  vortices 
was  taken  as  an  appropriate  correction  for  the  influence  of  the  portion  of 
wake  downstream  of  the  tenth  station  on  the  flow  upstream  of  it.  This 
method  was  employed  in  all  subsequent  stages. 

The  result  of  calculation  of  the  wake  shape  behind  a rectangular  wing 
of  AR  = 3 is  shown  in  Fig.  10.  A comparison  of  this  result  with  that  of  the 
infinite  vortex  (Westwater)  model  is  indicated  as  the  2-D  curve  in  this 
Figure . 

At  this  stage  a lifting-surface  theory  was  introduced  to  enable  us  to 
deal  with  general  wing  shapes  with  some  accuracy  (Fig.  4) . 

The  standard  vortex-lattice  technique  is  employed  [5  - 10]  to  start 
with.  The  wing  is  divided  into  cells  by  a number  of  equally  spaced  stream- 
wise  lines  and  by  either  equally-spaced  constant-percentage-chord  lines  or 
by  lines  perpendicular  to  the  stream.  A bound  vortex  element  is  placed  at 
the  "quarter-chord"  of  each  box  and  a control  point  is  established  at  the 
middle  of  its  "three-quarter-chord"  line.  The  usual  trailing  vortices  are 
taken  to  spring  from  each  bound  vortex  element  and  to  continue  straight 
downstream  in  the  wing  plane,  and  the  unknown  vortex  strengths  are  found, 
in  the  usual  way,  by  evaluating  the  downwash  at  each  point,  equating  it  to 
the  value  given  by  the  boundary  condition,  and  solving  the  resultant  system 
of  linear  equations.  The  wing  aerodynamic  coefficients  are  then  found  from 
the  results,  to  a linear  approximation. 

A first  approximation  to  the  rolling  up  is  now  found  by  taking  the 
actual  spanwise  distribution  of  the  circulation  (which  will  generally  not  be 
elliptic  and  performing  a two-dimensional  calculation  of  rolling  up,  as  was 
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done  in  Section  (a)  for  the  elliptic  case.  This  now  taken  as  the  first 
approximation  to  the  rolled  up  wake  shape  starting  from  the  trailing 
edge.  As  in  Section  (b) , the  first  approximation  to  the  wake  shape  is 
improved  iteratively  using  the  calculated  vortex  strengths.  During  these 
iterations  the  vortex  strengths  on  the  wing  are  kept  constant.  Using  the 
new  wake  shape  the  vortex  distribution  on  the  wing  and  the  aerodynamic 
coefficients  are  recalculated.  We  then  return  to  the  reclaculation  of  a 
new  vortex  wake  shape  by  iterations;  keeping  the  vortex  strengths  fixed, 
and  so  on.  Some  wake  shapes  calculated  by  these  procedures  are  shown 
in  Figs.  11  and  12. 

In  the  program  of  investigation  of  the  wing  with  vortices  leaving  the 
planform,  the  calculation  starts  with  the  free  vortices  leaving  at  an  angle 
of  a/2  to  the  wing  planform  and  then  extending  in  straight  lines  into  the 
wake . The  system  of  bound-vortex  elements  is  set  up  exactly  as  in  the 
previous  section  and  the  vortex  strengths  and  wing  properties  are  found  by 
a similar  process.  The  rolling  up  of  the  vortex  wake  can  now  be  obtained  by 
calculating  the  induced  motion  of  the  individual  vortex  lines  due  to 
mutual  interactions.  The  final  shape  is  determined  by  an  interation 
procedure  described  in  Fig.  5,  involving  the  recalculation  of  the  vortex 
strengths  and  the  corresponding  wake  shape.  This  calculation  cannot  be 
done  using  the  ideal  line  vortices  because  of  the  difficulties  of  inter- 
section and  "escape"  of  the  vortices  discussed  previously.  Therefore, 
this  model  is  used  only  to  determine  the  non-linear  aerodynamic  coefficients 
with  straight  trailing  vortices.  The  use  of  finite  core  vortices  may 
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eliminate  some  of  the  problems  and  enable  the  calculation  of  this  rolled- 
up  vortex  wake  to  be  carried  out.  This  calculation  is  now  underway.  How- 
ever a simpler  calculation  is  possible  when  it  is  assumed  that  vortices 
are  shed  only  from  the  planform  edges.  In  this  case  the  new  vortex  wake 
shape  can  be  recalculated  by  the  iterative  procedure  indicated  in  Fig.  6 
step  1 . 

The  numerical  problem  associated  with  the  use  of  ideal  line  vortices 
lead  us  to  look  for  a solution  by  the  use  of  finite  core  vortices.  Since 
the  actual  wake  has  thickness  and  is  not  an  infinitesimally  thin  sheet, 
this  is  also  a better  approximation  to  the  real  wake.  We  assume  a cere 
diameter  which  is  determined  by  the  condition  of  constant  vorticity  per 
unit  cross-sectional  area  of  the  core.  The  density  of  the  vorticity  m 
the  core  can  then  be  calculated,  for  instance,  by  using  the  core  size 
calculated  by  Spreiter  and  Sacks  [19] . 

We  impose  the  condition  that  as  soon  as  two  cores  touch  the  two 
vortices  are  replaced  by  a single  one  at  the  center  of  gravity  of  the  pair . 
The  new  core  diameter  is  determined  by  the  condition  of  conservation  of 
vorticity  at  the  same  votex  density.  This  procedure  overcomes  the  difficul- 
ties caused  by  the  approach  of  ideal  vortices,  i.e.,  "intersection"  and 
"escape".  This  procedure  limits  the  number  of  vortices  automatically, 
thus  helping  to  eliminate  also  other  irregularities.  Furthermore,  this 
method  ensures  that,  when  all  the  vortices  are  amalgamated  far  downstream, 
we  end  up  with  a pair  of  vortices  with  the  correct  diameter  and  strength 
matching  the  expected  vortex  pair  in  the  far  field. 
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Under  these  assumptions  the  core  radius  is  determined  by  the 
relation. 


R2 


1 


Z K. 

j-i  : 


o 


(2) 


Rq  being  the  final  core  diameter. 


The  vortex  core  diameter  is  evaluated  at  present  by  using  the 
method  presented  by  Spreiter  and  Sacks  [19]  . Accordingly  the  core 
diameter  is 

4tt  • AR  • CD 

= 2{1  + exp  [ ^ ” ll}-1  (3) 

CL 

2 

when  the  wing  loading  is  elliptical,  C = C AAR  and  b'/b  = ff/4. 

L) . L 

i 

In  this  case  the  core  radius  of  the  vortex  in  the  far  field  is 


R = 0.197b'  = 0.155b. 
o 

The  vortex-wake  shape  obtained  when  41  finite-core  vortices  are 
initially  placed  with  equal  spacing  on  the  trailing  edge  of  a rec- 
tangular wing  is  shown  in  Fig.  14.  It  should  be  pointed  out  that 
during  the  wake  development  calculations  a number  of  vortices  combine 
and  we  find  less  vortices  in  the  wake  than  the  initial  number. 


The  procedure  for  the  calculation  of  the  wake  characteristics 
associated  with  the  finite-core  vortices  is  outlined  in  Fig.  6.  The 
detailed  procedure  for  the  calculation  of  the  vortex  core  diameter  is 
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oulined  in  Fig.  7 
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3.1.  The  Thick  Vortex  Layer  Model 

This  investigation  is  an  attempt  to  remove  the  unrealistic  features 
of  the  earlier  models,  namely  zero  wake  thickness  and  vorticity  con- 
centrated on  lines,  by  assuming  that  the  wake  vorticity  is  contained  in  a 
layer  of  finite  thickness  with  some  plausible  cross-sectional  shape . The 
wake  flow  is  assumed  to  be  two-dimensional  and  the  rolling-up  is  studied 
via  the  time-dependent  development  of  this  model,  exactly  as  in  Westwater's 
work,  so  that  we  deal  with  a slow  rolling-up  taking  place  far  behind  the 
wing,  once  more.  The  introduction  of  vorticity  distributed  continuously 
throughout  the  wake  cross-section  enables  us  to  obtain  a more  acceptable 
picture  of  the  initial  rolling  up  phase  than  the  infinite  spiral  of  Kaden. 

Having  selected  an  appropriate  wake  cross-section  its  area  is  divided 
into  triangular  elements  within  each  of  which  the  vorticity  is  assumed 
constant,  and  for  which  simple  expressions  giving  the  velocity  field  have 
been  derived  (Ref.  10) . The  strength  of  the  vorticity  within  each  triangle 
is  determined  using  two  assumptions: 

(1)  The  vorticity  is  constant  trough  the  wake  thickness.  This 
corresponds  to  an  assumption  that  the  transverse  velocity 
profile  within  the  wake  is  linear. 

(2)  The  periphery  of  the  wake  is  moving  downwards  with  a velocity 
determined  by  spanwise  position  and  wing  spanwise  loading, 
exactly  as  in  ordinary  wing-wake  theory. 

Assumption  (1)  can  evidently  be  removed  at  the  cost  of  increasing 
the  number  of  triangular  elements  used. 
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Once  the  triangle  strengths  are  found,  the  network  of  points  defining 
the  wake  is  allowed  to  distort  with  time  under  its  self-induced  velocity 
field,  using  Euler  integration.  During  this  process  viscous  dissipation 
is  neglected,  so  that  the  vorticity  inside  each  triangle  remains  constant, 
as  will  its  area  (due  to  continuity),  even  though  the  shape  changes. 

The  mathematical  model  of  the  wake  and  the  notation  are  illustrated 
in  Fig.  15. 

The  unrolled  wake  cross-section  A' , whose  boundary  is  denoted  by 
C',  stretches  between  y = -s  and  y = +s  in  the  plane  x = 0.  At 
subsequent  stations  x = x,  the  boundary  is  denoted  by  C and  the  wake 
has  a rolled  up  form  of  cross  Section  A with  vorticity  distribution 
£(x,y,z)  determined  by  the  original  configuration  of  the  wake  cross- 
section,  the  original  vorticity  distribution  with  it,  £(0,y,z),  and  the 
elapsed  time  x/U.  It  is  assumed  that 

£(0,y,z)  = co  (y)  , (4) 

corresponding  to  a linear  variation  with  z of  the  velocity,  v,  within 
the  initial  wake  section,  if  the  wake  is  assumed  thin. 

We  define  dimensionless  coordinates 

* 


and  since  a typical  velocity  of  the  flow  in  the  wake  cross-sectional 

plane  is  2UC  /w *AR,  we  also  define  dimensionless  velocity  components, 
L 


time  and  vorticity  by 
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* ir  AR.v 

v = 

2UC 

L 

t*  2Clx 

TT  . AR . S 


* ttAR.W 
W " 2UC 

Li 


*,  * ojtt.AR.s  * * * * Stt.AR.s 

“ (y  > = IOC J 5 (t  'y  '*  )=  ~ 2UC~ 

L L 


We  assume  that  oo  or  w is  determined  by  the  boundary  con- 


ditions : 


on  C’ 


where 


* it  * * 

w cosX  - v sinX  = w (y  ) cosX 


* * 

(dz  /dy  ) C * = tanX,  and 


*2  * 
v + w 

\ 


O as  y + z 


- w (y  ) is  the  non-dimensionalised  downwash  distribution  in  the 

Trefftz  plane  as  calculated  from  the  usual  unrolled  thin-wake  theory. 

* 

For  example,  = - 1 for  elliptic  spanwise  load  distribution.  It 

* 

is  readily  shown  that,  in  the  general  form,  w^  is  dependent  on  the 

★ 

form  of  the  spanwise  loading  and  y only  (c.f . Reference  3) . 

It  now  follows  that  for  a given  initial  wake  cross-section  and  span- 

★ * 

wise  loading,  the  velocities  v and  w of  a given  fluid  particle 


are  functions  of  t only,  so  that  the  subsequent  non-dimensional  co- 
ordinates of  the  particles  constituting  the  rolled  wake  are  functions  of 


t only,  determined  by  the  differential  equations 


The  equations  (8)  are  integrated  numerically  step  by  step,  starting  from 


the  initial  configuration,  using  Euler  integration. 
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It  remains  to  determine  the  velocity  field  (v  , w ) due  to  the  distribution 
* 

£ within  C. 


At  any  value  of  x,  the  two-dimensional  velocity  field,  (v,  w) , is 
determined  bv  integrating  the  effects  of  point  vortices  of  strengths 
£(x,y,z)  dy  dz  over  the  area  of  the  wake  cross-section,  A. 


To  facilitate  the  numerical  calculations,  A is  divided  into  a 

finite  number  of  small  triangular  elements  within  each  of  which  the  value 

* ★ 

of  £ is  assumed  constant.  We  now  find  expressions  for  Av^  and  Aw^ 
due  to  a typical  element  of  this  kind  for  points  outside  its  boundary,  or 
approaching  the  boundary  in  a limiting  sense.  This  enables  us  to  calculate 

* it 

v and  w for  the  assembla^  of  triangles  at  all  the  node  points  of  the 
tringular  mesh,  including  internal  nodes  of  A,  since  such  points  may  be 
regarded  as  being  inside  infinitesimal  cavities  excluded  from  all  the 
adjacent  triangles,  and  we  calculate,  in  effect,  the  principal  value  of 
the  velocity  integral  - which  is  precisely  the  required  definition  of  this 
integral  inside  the  vorticity  distribution. 


The  velocity  field  outside  any  area  S containing  a constant  vorticity 

distribution,  £ , may  be  written 
s 


v 


iw  = 


-’s  dy^dz^ 

-s  y + iz  - - izi 


(12) 


Using  Green's  theorem,  this  may  be  converted  to  a line  integral  around  P, 
the  perimeter  of  S: 


r 
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v - iw  = 


2tt 


£n(c  ~ Y1  ~ 


(13) 


where  £ = y + iz. 

If  (10)  is  applied  to  a typical  triangle  of  A,  of  area  AA  , 


vorticity  distribution  strength  £ and  vortices  defined  by  c ,r 

n 1 A 


Cc  (Fig.  2) , we  obtain 


£ - £, 


-i£  A A 

Av  - iAw  = 2 — H {_ - 

n n * (£-£*)(£-£,) 


C A'  B A' 

£ - £, 


*n(£  - Cft)  + 


S -cB  _ , ,c 

<£„-£„)  <£.-£„>  inU~iB)  + (£  -c  )U  -c  ) 


C B A B 


A ’C  B ’C 


An(C-Cc)} . 


(14) 


£ and  Aa  are  both  invariant  with  the  motion  for  an  infinitesimal 
n n 


triangle  in 

incompressible  flow. 

the  product  £ AA 
n i 

culation  around  the 

element . 

If  we  now  write 

2 * 

s Aa  = aa 

, then 

£ AA  = 

n n 

2UsCL  * * 

— ■ £ AA 

2 UsC 

= — - K 

n n 

TTAR  n n 

irAR  n 

and  if  (see 

Fig.  2) 

£ - £, 

ia 

= sr  e ; 

£ - £ = sr  e^ 

A 

A 

B B 

S^A 

iaB 

= see  D ; 

‘c  - ■ s-iaC 

i$,-i 

= sae  ; 

_ r i6A 

£A  - £B  = see 

= sbe^Y A 

r r iyB 

£b  - £c  = sae 

(15) 


£ - £c  = srce 


1Y 


(16) 


we  obtain  from  (11) 
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* Kn  rA 

AVn  = 71  {b^[sin(0t  ~ “c  " aBUnrA  + aCOS  (0t  " aC  “ aB)J  + 


1 


+ — [since  - ec  - 6A)HnrB  + Bcos(B  - Bc  - Bft) ] + 


+ Sh[sin(Y  ~ ya  “ YBUnrc  + ycos(y  “ Y*  - Yj]}  = 


= Vn^'2*' 


A ’B 


(17) 


* Kn  rA 

Awn  = — ^ [cos (a  - ac  - oig)  £nrA  - asin(a  -ac  - ag)  ] 


+ — [cos(B  - 6C  - 6A)inrB  - Bsin(B  - BQ  - BA> ] + 


+ ab  [cos(y  " Ya  * Yn)lnrr  ~ Ysin(y  - Y*  - Y J ] > = 


’A  'B'  C 


A 'B 


= K W (y*,z*) 

n n 


(18) 


Hence  we  have  the  approximations 

M ^ „ 


% 


Z K V (y  ,z  ) 
, n n 
n=l 


(19) 


w % Z K W (y  ,2  ) . 
, n n 
n=l 


(20) 


The  details  of  the  numerical  calculation  procedure  and  results  for 
the  wake  of  elliptic  cross-section  are  presented  in  Ref.  10.  The  notation 
for  the  triangular  elements  is  shown  in  Fig.  16.  Some  results  showing 
the  appearance  of  the  tip  region  for  various  stages  of  the  initial  roll- 
up are  shown  in  Figs.  18,  19  and  20  for  e = 0.4,  0.5  and  0.6,  respectively 


— - — - — — 


i 


and  in  Fig.  21  one  example  of  a complete  half-wake  is  shown.  It  should  be 
noted  that  to  clarify  the  inner  detail  of  the  spiral,  the  vertical  scale 
has  been  exaggerated  in  these  figures  and  this  vertical  scale  is  not 
uniform  throughout  all  the  pictures.  The  results  of  this  calculation  and 
the  significant  features  of  the  model  used  are  also  discussed  in  detail 
in  Ref.  10. 

The  analysis  of  the  results  shows  that,  for  the  range  of  thickness 
ratios  covered,  for  the  elliptic  cross-sectioned  wake: 

(1)  Thickness  has  little  effect  on  the  amount  of  vorticity  within 

★ 

the  core  at  a given  t (Figure  22) . 

(2)  Thickness  has  almost  no  effect  on  the  size  of  the  core  at  a 

* 

given  t (Figure  23) . 

(3)  Except  in  the  very  initial  stages  of  roll-up,  the  core  is  of 
approximately  elliptic  shape  with  height/width  ratio  about  0.8 
for  all  three  thickness  ratios  studied  (Figure  24) . 

(4)  The  number  of  turns  of  the  spiral  within  the  core  at  a given 

* 

t increases  with  reduction  of  e,  except  for  very  small 
* 

values  of  t when  the  spiral  is,  in  any  case,  ill-defined 

(Figure  25) . About  one  and  three-quarter  turns  are  observed 

* 

for  the  e = .04  case  at  t = .0128  and  for  the  e = .05  case 

* * 
at  t = .0160,  and  for  the  e = .06  case  at  t = .0192,  just 

over  one  and  a half  turns. 

(5)  An  interesting  result  observed  in  all  three  cases  is  that  the 
original  tip  point  of  the  unrolled  wake  does  not  become  the 


24 


i 


tip  point  of  the  spiral;  instead  it  recedes  back  along  the 
outer  edge  of  the  coil.  This  appears  to  be  true  also  for 
subsequent  tip  points  of  the  spiral  - these  do  not  remain  at 
the  tip  but  are  dragged  back  along  the  outer  edge  of  the 
spiral  in  their  turn  and  their  place  is  taken  by  other  points 
which  were  originally  more  inboard  along  the  upper  edge  of 
the  unrolled  wake . 

This  work  is  being  extended  into  the  more  complete  stages  of  rolling- 
up  and  also  to  deal  with  different  spanwise  loadings. 


I 
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IV.  THE  METHOD  OF  CALCULATION  OF  THE  NONLINEAR  AERODYNAMIC  CHARACTERISTICS 
OF  WINGS. 

] 

The  present  method  is  developed  to  calculate  the  aerodynamic  charac- 
teristics of  wings  at  subsonic  speeds.  The  wings  can  have  various  Aspect 
Ratios  from  low  to  high  and  can  be  at  any  angle  of  attack  below  that  of 
stall  or  vortex  breakdown.  The  calculation  includes  also  the  evaluation 
of  the  shape  of  the  rolled  up  vortex  sheet  above  and  behind  the  wing  as 
well  as  the  velocity  field.  A schematic  model  for  the  flow  over  a wing  at 
a high  angle  of  attack  as  presented  by  Maltby  [34]  is  shown  in  Fig.  26. 

For  low  aspect  ratio  wings,  the  rolling  up  of  the  vortex  sheet  above  the 
wing  affects  the  pressure  distribution  and  results  in  the  non-linear  lift 
variation  at  increasing  angles  of  attack. 

The  calculation  procedure  is  described  in  the  flow  chart  shown  in 
Fig.  1.  The  calculation  is  divided  into  the  following: 

a.  Definition  of  the  Model  for  the  Wing  Planform  (Steps  1 and  2) 

The  planform  of  the  wing  is  divided  into  cells  as  shown  in  Fig. 

2.  The  number  of  subdivisions  has  to  be  limited  due  to  nume- 
rical problems  as  well  as  the  size  of  computer  memory.  The  bound 
vortex  and  the  control  points  positions  are  clearly  indicated  in 
Fig.  2. 

b.  Establishment  of  the  Initial  Conditions  and  the  Strength  of 
the  Vortices  (Step  3)  . 

The  initial  conditions  for  the  calculation  involve  the  deter- 
mination of  the  values  of  the  strength  of  the  vortices  in  each 


cell  on  the  planform  and  in  the  wake.  This  requires  the 


™r~ 
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assertion  of  an  initial  shape  for  the  wake.  For  rectangular 
and  delta  wings  the  initial  shape  of  the  vortices  is  assumed 
to  be  straight  lines  issuing  from  the  bound  vortex  of  each 
cell  and  continuining  downstream  towards  + » . These 
vortices  can  be  issued  from  the  trailing  edge  only, so  that 
initially  the  vortices  are  embedded  in  the  planform.  Other 
options  are  to  issue  the  vortices  in  straight  lines  at  an 
arbitrary  local  angle  from  : (a)  each  cell,  (b)  from  all  edges 
(sides,  leading  and  trailing  edges)  or  (c)  from  side  and 
trailing  edges  only.  These  options  can  be  selected  by  the 
proper  computer  code.  The  initial  local  angle  is  selected 
to  be  half  the  angle  of  attack.  This  trailing  vortex  con- 
figuration which  was  suggested  by  Bollay  [35]  and  by  Gersten 
[36]  is  found  to  be  very  suitable  as  a starting  condition. 

These  initial  conditions  enable  the  calculation  of  the  correspond- 
ing strengths  of  the  vortices. 

c . Calculation  of  the  Velocity  Field  (Step  4) . 

The  induced  velocity  at  each  point  is  calculated  by  the  Biot- 
Savart  law,  by  adding  the  influence  of  all  the  vortices  as  well 
as  self  induced  effects.  The  effect  of  curvature  on  the  self 
induced  velocity  is  dealt  by  choosing  a cutoff  radius  as  suggested 


by  Crow  (Ref.  32) . 


* 

* 
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d.  Calculations  of  the  Vortex  Trajectories  (Step  5) 

The  position  of  the  vortices  is  calculated  from  the  induced  flow  field 
by  the  numerical  integration  of  the  relations  presented  in  Part  II,  Eq. 

1.  During  the  investigation  both  the  Euler  integration  method  and  a 
second  order  Runge-Kutta  method  were  used. 

In  order  to  facilitate  the  convergence  of  the  calculation  the  iteration 
procedure  for  the  determination  of  the  new  vortex  position  is  damped  by 
a relaxation  factor. 

e.  Calculation  of  the  Strength  of  the  Vortices  (Step  6) 

Knowing  the  shape  and  position  of  the  vortices  in  a given  iteration  the 
strengths  of  the  vortices  are  calculated  by  a solution  of  the  matrix 
equation 

{tan  a}  = [H]  {k}  (21) 

| 

where  [H]  is  the  influence  matrix  and  (k)  is  the  strength  of  the 
vortices  while  a is  the  local  angle  of  attack.  The  order  of  the 
matrix  is  equal  to  the  number  of  cells.  Each  equation  is  obtained 
by  imposing  the  zero  vertical  velocity  at  the  control  point  of  each 
cell.  The  matrix  equation  (21)  is  solved  numerically. 

f . Calculation  of  the  Aerodynamic  Coefficients  and  the  Wake  Characteristics 
(Step  7) 

At  the  completion  of  each  iteration  the  vortex  strength  and  the 
positions  of  the  vortices  are  determined.  These  vortex  calculations 
can  be  used  for  the  evaluation  of  the  aerodynamic  loading  distribu- 


tions as  well  as  the  integral  values  for  lift,  drag  and  various 


It  is  interesting  to  note  in  Fig.  27  that  the  main  features  of  the 
vortex  sheet  shape  are  already  obtained  when  the  vortices  are  allowed  to 
leave  the  wing  at  the  side  edges  only.  Furthermore,  there  is  practically 
no  difference  between  the  results  obtained  when  vortices  were  allowed  to 
leave  from  all  cells  or  from  only  the  leading  and  side  edges. 

The  interesting  effects  of  non-linear  loading  variation  can  be  seen 
from  the  calculations  of  the  loading  on  a rectangular  wing  of  AR  = 1 . 

This  wing  is  selected  because  of  the  available  experimental  data  reported 
in  Refs.  39,  40  and  41. 

The  results  of  the  calculations  using  the  present  method  and  the  com- 
parison with  the  experimental  data  for  C^,  and  CD  are  shown 

i 

in  Figs.  29,  30  and  31  respectively.  The  planform  divisions  used  are: 

N = 2,  N =10;  N = 2,  N =20;  N = 4,  N = 10.  For  all  of  these 
c s c s c s 

divisions  the  results  are  very  good.  The  value  of  CQ  increases 

i 

slightly  for  larger  Ng,  improving  the  agreement  with  the  experimental 

data.  In  Fig.  31,  the  calculated  values  of  CD  are  shifted  by  the 

i 

experimentally  measured  C for  better  comparison  with  the  measured 

D0 

values . 

Comparison  of  the  experimental  and  calculated  loading  at  a = 20° 


and  a = 10  are  presented  in  Figs.  32,  33,  34  and  35. 


Here,  the  agreement  is  generally  good  except  for  the  high  peaks  in 
ACp  near  leading  and  side  edges.  This  is  well  illustrated  in  Fig.  36 
where  the  effect  of  the  wings'  profile  measured  in  Ref.  40  are  presented. 

It  can  be  seen  that  the  calculated  results  are  in  better  agreement  with 
the  data  of  the  thin  wing. 

The  shape  of  the  trailing  vortex  wake  is  shown  in  Figs.  37,  38  and 
39.  The  roll  up  of  the  vortices  at  the  trailing  edge  plane  of  the  wing 
is  illustrated  in  Fig.  37.  It  can  be  seen  that  the  vortex  trajectories 
are  such  that  there  is  a large  vortex  on  the  wing  situated  near  each  of 
the  side  edges  over  the  planform.  The  interaction  of  the  rolled  up  side 
vortices  with  the  vortices  shed  from  the  trailing  edge  is  shown  in  Fig.  38. 
The  rolled  up  vortex  contains  42%  of  the  total  wing  vortex  strength.  At 
the  distance  of  0.8b  behind  the  wing  this  vortex  is  much  larger  in  size 
and  in  strength.  The  vortices  shed  from  the  trailing  edge  are  much 
weaker.  The  vortices  near  the  edge  are  of  opposite  sign  to  those  at  the 
rolled  up  side  vortex.  Therefore,  a cusp  is  obtained  at  their  inter- 
section. The  observation  of  Maltby  (Ref.  34),  shown  in  Fig.  26,  indeed 
shows  the  appearance  of  this  cusp  and  the  existence  of  a secondary 
vortex  of  opposite  sign  at  this  intersection  position.  The  overall  shape 
of  the  vortex  over  and  behind  the  AR  = 1 wing  up  to  10  half  spans  behind 
the  trailing  edge  is  shown  in  Fig.  39. 

The  availability  of  experimental  results  for  a rectangular  wing  of 

AR  = 0.25  (Ref.  42)  make  this  wing  a good  test  for  our  method  for  the  very 

slender  rectangular  wings.  The  aerodynamic  coefficients  are  compared  in 

o 

Excellent  agreement  is  found  even  at  a = 20  . 


Fig.  40. 
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The  pressure  loading  , tc.^  , on  the  wing  is  shown  in  Figs.  41  and  42, 
the  measured  and  calculated  values  are  presented  respectively.  The  agree- 
ment between  these  results  is  very  good  except  at  the  leading  edges. 

These  calculations  are  performed  with  the  wing  divided  into  48  cells 

(N  = 8,  N = 6) . This  division  was  selected  in  order  to  facilitate  the 
c s 

comparison  with  the  experimental  data.  The  shape  of  the  rolled  up  vortex 
can  be  obtained  from  the  computed  trajectories  of  the  side  vortices  shown 
in  Fig.  43. 

(2)  DELTA  WINGS 

The  triangular  geometry  of  the  delta  planform  requires  special  con- 
siderations due  to  the  cell  shape  and  therefore  to  the  applicable  cal- 
culation method. 

There  are  a number  of  methods  proposed  in  Refs.  17,  22  and  23,  how- 
ever, we  are  able  to  obtain  good  results  with  the  computation  method 
presented  in  Figs.l  and  2.  The  triangular  planform  is  divided  into  equal 

number  of  chorwise  and  spanwise  divisions,  N = N . The  triangular  cells 

c s 

have  their  bound  vortex  on  a line  connecting  the  l/4c  point  to  the  tip. 

The  control  point  is  located  on  the  median  line  at  the  3/4c  distance  from 
the  apex . 
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The  calculated  lift  coefficients  for  delta  wings  of  AR  = 1;  1.5  and  2 
are  compared  with  the  experimental  data  as  well  as  with  the  classical  non- 
linear theories.  These  are  shown  in  Figs.  44,45  and  46  respectively  The  agree- 
ment with  the  experimental  results  is  reasonable,  however,  the  variance 
increases  as  the  aspect  ratio  decreases.  The  aerodynamic  center  for  the 


AR  = 1 delta  wing  is  shown  in  Fig.  47. 


The  loading  in  each  cell,  AC^,  is  calculated  and  compared  to  the 
experimental  data  of  Ref.  47  in  Figs.  48  and  49  at  nominally  a = 20° 
and  at  a = 14°  in  Figs . 50  and  51 . The  calculated  distributions  shown 
in  Figs.  49  and  51  are  similar  to  the  measured  ones  although,  there  is  loss 
of  details  in  the  distribution  due  to  the  relatively  large  cell  size  and 
due  to  the  edges  effects  noted  already  in  the  treatment  of  the  rectangular 
wings . 

The  induced  drag  is  also  evaluated  and  compared  to  the  experimental 
results  of  Ref.  48  in  Fig.  52.  Here  the  agreement  is  very  good. 


The  pattern  of  the  vortices  shed  from  the  delta's  leading  edges  is 

presented  in  Fig.  53  for  a wing  of  AR  = 1 and  a 6 x 6 division  (N  =N  =6) . 

c s 

Although,  the  results  presented  here  are  encouraging  there  are 
difficulties  in  the  calculations  of  delta  planforms,  particularly  for 
slender  wings  with  AR  less  than  1.  These  difficulties  arise  because  of 
the  very  strong  induced  velocities  due  to  the  small  distances  between  the 
discrete  vortices  at  low  AR  and  large  number  of  cells.  If  one  introduces 
finite  core  vortices,  the  core  size  limits  the  number  cells.  Therefore 
it  is  reasonable  to  expect  good  results  with  the  present  method  for 
delta  wings  of  Aspect  Ratios  1,  or  higher. 

(3)  SECONDARY  VORTEX  SEPARATION 

The  present  method  includes  the  option  of  calculating  the  trajectories 
of  the  vortices  when  these  are  shed  from  each  cell.  This  calculation 
results  in  the  appearance  of  additional  rolled  up  vortices.  These  secondary 
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vortices  are  much  weaker  than  those  that  form  the  main  lift  vortex.  The 
existence  of  these  secondary  vortices  is  found  in  experiments  on  both 
rectangular  wings  and  on  delta  wings  (Ref.  34) . The  calculations  of  the 
secondary  vortices  on  delta  wings  involved  some  problems  due  to  one  of  the 

i 

"secondary"  vortices  jumping  outside  of  the  wing's  planform.  However,  here 
also  the  secondary  vortex  phenomenon  is  observed. 

(4)  COMPULATION  PERIODS 

The  present  method  requires  about  150-800  seconds  on  the  IBM  370/168 
computer . 

The  rectangular  wings  calculations  require  about  150-200  seconds, 
while  the  delta  wing  calculation  require  more  iterations  so  that  the  com- 
puting time  is  extended  to  about  300  seconds.  The  calculations  of  the 
full  vortex  wake  trajectories  more  than  doubles  these  periods,  to  about 
600-800  seconds . 
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V.  DISCUSSION  OF  RESULTS 

The  present  method  is  based  on  the  extension  of  the  Vortex  Lattice 
concepts  also  to  low  aspect  ratio  wings  and  to  the  vortex  sheet  which  is 
established  over  and  behind  all  wings.  As  part  of  this  investigation  a 
computer  program  is  developed.  This  program  is  capable  of  reproducing 
the  results  of  the  Vortex  Lattice  Method  (VLM)  for  linear  aerodynamic  co- 
efficients of  high  aspect  ratio  wings.  The  addition  of  the  effect  of  the 
vortex  sheet  behind  the  wing  (i.e.  vortices  shed  from  the  trailing  edges) 
is  presented  in  the  Modified  Vortex  Lattice  Method  (MVLM) . The  program 
enables  the  calculation  of  the  rolled  up  vortex  sheet  behind  the  high 
aspect  ratio  wings  as  well  as  the  correct  lift  (circulation)  distribution 
on  the  wing  planform.  This  wake  calculation  showed  that  for  rectangular 
wings  the  center  of  the  tip  vortex  core  is  very  close  to  the  wings ' tip 
rather  that  at  the  n/4(b/2)  obtained  from  the  elliptic  lift  distribution. 
Furthermore,  it  was  found  that  although  the  vortex  wake  behind  the  wing 
rolls  up  into  two  concentrated  vortices  practically  these  do  not  affect 
the  linear  lift  of  the  wing  itself.  It  is  therefore  clear  that  the  non- 
linear aerodynamic  coefficients  are  due  to  the  vortices  which  are  shed  on 
and  over  the  wings'  planform.  This  effect  is  particularly  significant  for 
low  aspect  ratio  wings.  Therefore  the  calculation  method  does  include 
the  programs  for  the  evaluation  of  trajectories  of  the  vortices  shed  from 
the  side  and  leading  edges  and  from  each  cell  on  the  planform.  These 
vortices  roll  up  into  the  well  known  leading  edge  vortices  on  the  delta 
wing  or  to  Die  attached  side  vortex  for  the  low  aspect  ratio  rectangular 
wing.  The  resulting  non-linear  aerodynamic  characteristics  are  computed 


with  good  accuracy. 
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The  representation  of  the  vortex  sheet  by  discreet  vortices  causes 
problems  in  the  numerical  calculations.  It  was  found  that  the  number  of 
divisions  has  to  be  limited.  In  order  to  allow  for  a more  realistic  model 
the  finite  core  vortex  model  has  been  proposed.  Some  calculations  using 
this  model  have  also  been  performed. 

The  results  of  the  method  in  the  evaluation  of  the  aerodynamic  charac- 
teristics on  wing  and  their  wakes  are  very  encouraging.  By  allowing  vortices 
to  be  shed  from  all  the  cells  in  the  planform  the  vortex  structures  that 
are  obtained  in  these  calculations? the  secondary  vortex  separations,  are 
similar  to  those  found  in  experimental  measurements. 

The  calculations  of  the  pressure  distributions  for  both  rectangular 
and  delta  wings  are  generally  in  good  agreement  with  experimental  data. 

Some  differences  are  due  to  the  local  effects  of  the  leading  edge  shape 
which  affects  the  results  of  the  experiments  • These  leading  edges  shape 
effects  are  not  accounted  for  in  the  present  calculation  method. 

The  extremely  short  distances  between  the  vortices  over  the  low  aspect 
ratio  delta  wings  causes  considerable  problems,  it  was  found  that  the  present 
method  can  be  used  for  the  computation  of  delta  wings  down  to  AR  = 1.  Some 
difficulties  appear  with  the  secondary  vortex  structure,  where  some  of 
the  vortices  are  "thrown  out"  of  the  vortex  pattern  around  the  wing.  These 
effects  must  be  further  studied. 

In  conclusion,  it  seems  that  the  present  method  is  capable  in  dealing 


with  a large  group  of  wing  planforms  and  present  results  which  are  in  good 
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agreement  with  experimental  data.  It  is  planned  to  extend  the  computer 
programs  for  planforms  of  general  shapes  and  to  include  effects  of  camber 
and  flaps . The  method  can  also  be  extended  to  the  calculation  of  several 
wings  in  tandem,  thus  enabling  the  evaluation  of  the  aerodynamic  charac- 
teristics of  the  mutually  interactive  wings. 
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FIG.  1 -FLOW  CHART  OF  CALCULATION  METHOD 


FIG.  2 -THE  CELL  ARRANGEMENTS  AND  THE  ASSOCIATED 

VORTICES  AND  CONTROL  POINTS  FOR  RECTAN- 
GULAR AND  DELTA  WINGS. 
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FIGURE  3 - VORTEX  WAKE  CALCULATIONS  FOR  ELLIPTIC  LIFT  DISTRIBUTION 
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FIGURE  13  - THE  ROLLED- UP  WAKE  SHAPES  CALCULATED 
STEP  SIZE  M = 16  AND  M = 32. 
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FIG.28- COMPARISON  OF  THE  MEASURED  CORE 
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FIG.  31 -THE  INDUCED  DRAG  COEFFICIENT  OF  A RECTAN 
GULAR  WING  OF  AR=1. 


FIG.  32- ACp  OF  A RECTANGULAR  WING  OF  AR  = 1 AT  oc  = 20°  • EXPERI- 
MENTAL RESULTS  OF  SCHOLZ  [40]. 


ACp  OF  A RECTANGULAR  WING  OF  AR=1  AT  <*=10°,  EXPERI 
MENTAL  RESULTS  OF  SCHOLZ  [40]  . 


FIG.35-THE  ACp  CALCULATED  BY  THE  PROPOSED  METHOD  FOR 
A RECTANGULAR  WING  OF  AR=1  AT  or  = 10°. 
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FIG.38-THE  NEAR  WAKE  DEVELOPMENT 
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FIG.40  COMPARISON  OF  THE  AERODYNAMIC  COEFFI 
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WING  OF  AR  = 0.2  5 . 
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OF  AR=a25  AT  INCIDENCE  o(  = 20°. 
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F I G.44 -COMPARISON  OF  EXPERIMENTAL  AND  COMPUTED 
LIFT  COEFFICIENT  Ci  AS  A FUNCTION  OF  at 
FOR  DELTA  WING  OF  AR-1. 


O - EXPERIMENTAL  RESULTS  [A3] 


FIG. 45-COMPARISON  OF  EXPERIMENTAL  AND 

COMPUTED  LIFT  COEFFICIENT  Ci  AS  A 
FUNCTION  OF  o ( FOR  DELTA  WING  OF  AR-L5 


AR  = 2.0 


// 
i i 

// 


& $ 

f : 

y f 

/ ,0/ 

^ / 

o/  ^T/ 

// 


° ^ 

^ vl*/ 

J v" 


*vv 

VN> 

/ 


A - PRESENT  CALCULATION 
O- EXPERIMENTAL  RESULTS  [43j 


FI646  COMPARISON  OF  EXPERIMENTAL  AND  COMPUTED  LIFT 
COEFFICIENT  CL  AS  A FUNCTION  OF  cf  FOR  DELTA 
WING  OF  AR  s 2 . 
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SACKS  [46] 

O RE H BACH  [17] 

□ PECKHAM  EXPERIMENTAL  RESULTS  [45] 
• PRESENT  WORK 
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FIG. 47 -CENTER  OF  PRESSURE  OF  A DELTA  WING 
OF  AR»  1 AS  A FUNCTION  OF  INCIDENCE. 


FIG. 49  - THE  ACp  CALCULATED  BY  THE  PROPOSED  METHOD 
FOR  A DELTA  WING  OF  AR  = 1.46  AT  «=20°. 


FIG.51  -THE  ACp  CALCULATED  BY  THE  PROPOSED  METHOD  FOR 
A DELTA  WING  OF  AR=1.46  AT  of=15°. 


FIG.52  LIFT  COEFFICIENf  Cf  VS  INDUCED  DRAG 
COEFFICIENT  CDj  (AR=1.46). 
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G.53-VORTEX  LINES  ON  A DELTA  WING  OF  AR»1  AT  of  =2  0° 
(6x6  SUBDIVISIONS  ). 


